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Abstract 

The  problem  of  computing  the  reflected  field  due  to  the  incidence 
of  a  plane  electromagnetic  wave  on  a  plane  stratified  inhomogeneous 
raediuin  is  considered*  An  approximate  solution  to  this  problem  is 
examined,  and  the  error  is  estimated.  It  is  found  that  when  the  raediuin 
is  slowly  varying,  in  a  sense  that  can  be  defined  precisely,  the  error 
is  small.  It  can  also  be  seen  that  the  error  will  be  small  in  this 
case  no  matter  what  the  wavelength  may  be. 
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lo   Introduction 

This  paper  is  concerned  with  the  problem  of  calculating  the  electro- 
magnetic reflection  properties  of  a  plane  stratified  inhomogeneous  dielectric 
medium  such  as  an  atmospheric  layer  of  the  type  indicated  by  certain  micro- 
wave refractometer  measurements.  We  shall  assume  that  the  electrical  pro- 
perties of  the  medium  do  not  vary  greatly,  which  is  apparently  the  case  in 
an  atmospheric  layer,  but  we  shall  allow  the  medium  to  have  a  variable 
finite  conductivity  since  no  essential  complication  occurs  when  we  do  so. 
The  relations  between  the  electrical  parameters  are  assamed  to  be  such  that 
the  phenomenon  of  total  reflection  cannot  take  place. 

An  approximate  method  for  calculating  the  reflected  wave  due  to  an 
incident  plane  wave  will  be  described,  and  the  error  committed  in  the  use 
of  the  approximation  will  be  estimated.  We  shall  show  that  this  approxi- 
mation is  correct  for  all  wavelengths  to  within  the  error  estimate  which 
is  independent  of  the  wavelength.  Moreover,  for  the  validity  of  the 
method,  the  variation  of  the  medium  is  not  required  to  be  smooth  but  only 
continuous  with  piecewise  continuous  derivatives. 

The  approximation  we  shall  use  has  been  given  in  one  form  or  another 
by  several  authors.  It  is  a  special  case  of  the  method  devised  by  J,  B, 
Keller  and  H.  B,  Keller'-  -•  for  a  system  of  differential  equations.  The 
same  approximation  was  given  by  Bremmer"-  -•  and  by  Schelkunof f  *■  -■ ,  and  it 
is  derived  in  the  thesis  of  Landau'-''-',  The  approximation  is  often  regarded 

r2i 

as  a  correction  to  classical  optics  in  one  dimension  or  the  W.K.B,  solution^  J, 
hence  the  implication  that  its  accuracy  is  frequency  dependent.  We  shall  see, 
however,  that  while  the  accuracy  may  increase  with  increasing  frequency  the 
error  can  be  bounded  by  a  small  number  which  is  independent  of  the  frequency. 

This  paper  is  divided  into  sections  so  that  a  part  or  all  of  the  analysis 
can  be  easily  skipped,  if  desired.  In  fact.  Section  7  contains  a  summary  of 
the  results  necessary  for  applications. 
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2.   Statenient  of  the  mathematical  problem 

The  description  of  oiir  problem  will  be  given  in  terms  of  a  Cartesian 
coordinate  system  which  has  its  2-axis  oriented  vertically  so  that  z  increases 
in  the  upward  direction.  The  reflecting  medium  is  assumed  to  be  stratified  in 
horizontal  layers,  all  electrical  properties  of  the  medium  being  constant  in 
any  plane  parallel  to  the  x,y-plane  of  the  coordinate  system. 

Since  the  electrical  properties  of  the  medium  can  vary  in  the  z  direction 
only,  the  dielectric  parameter  eCz),  the  permeability  |J.(z)  and  the  conducti\'ity 
<r(z)  can  be  func  ions  of  z  but  not  of  x  and  y.  It  is  assximed  that  the  variable 
part  of  the  medium  has  a  lower  boundary,  the  x,y -coordinate  plane.  Thus,  c,  n 
and  6  are  all  constant  for  z  <  0,  We  designate  the  electrical  parameters  for 
the  homogeneous  medium  in  the  region  z  <  0  by  e  ,  u  ,  (T  ,  The  inhomogeneous 
part  of  the  medium  either  extends  upward  to  infinity  or  else  it  is  boianded  by 
a  plane  z  «  A  some  positive  distance  [\^  above  its  lower  boundary  at  z  ■  0, 
For  z  >  ^  the  medium  is  assumed  to  be  homogeneous  again. 

A  plane  electromagnetic  wave  of  arbitrary  polarization  is  incident  at  an 
angle  a  from  the  homogeneous  region  z  <  0,  The  coordinate  system  if  fiirther 
oriented  so  that  the  y,z -coordinate  plane  is  the  plane  of  incidence. 

Maxwell's  equations  for  a  wave  of  circular  frequency  co  will  be  used  in 
M,K,S,  units: 

(2,])  y  X  F  -  ->mH 

V  X  H  «  ((T  +  Jwe)E  -  jco£  E, 

where  6  ■   (<5"  +   jojeVjco,   and  the  time  factor  has  been  tal<©n  to  be  exp(ja)t). 

It  is   convenient  to  think  of  the  electromagnetic  field  as  the  sum  of  two 
separate  fields:  one  with  the  E  vector  polarized  parallel  tc  the  plane  of 
incidence,  and  one  polarized  normal  to  the  plane  of  incidence.     Thus,  the 
field  represented  by  the  pair  E,  H  is  a  sum  of  the  fields  represented  by 
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the  pairs  E  ,  H  and  E  ,  H  ,  where  the  subscript  p  indicates  polarization 
parallel  to  the  plane  of  incidence  and  n  indicates  polarization  normal  to 
the  plane  of  incidence.  Since  the  y,z -coordinate  plane  is  the  plane  of 
incidence  and  a  is  the  angle  of  incidence  we  have: 

(2.2)     Ep  -  (0,Vp,6)exp(-Jk^y  sin  a),    Hp  -  (lp,0,0)exp(-jk^y  sin  a) 

^n  "  (^n'^'^^^'^^'^V  ^^"  °^'     "n  '  (0,I„,-^  )exp(-Jk^y  sin  a). 


where  k     ■  co  /c  \x      is  a  propagation  constant  referred  to  the  homogeneous 

0        CO 

medium  of  the  region  z  <  0. 

From  these  assumed  forms  for  E  »H  ,F  ,H  ,  by  substituting  into  (2.1),  we 
obtain  the  two  pairs  of  transmission  line  equations: 

(2.?)         dV^/dz     -     -JoiL^Ij^ 

dl^/dZ     -     -Jo^C^V^,  i  -  P,n 

?  A  2 

with  L     -   (li  c  sin  a)/  ^  -  n,  C     -  -c,  L     ■  |i,  C     ■  e  -  (u  e  sin  a)/|J..     We 
also  have  the  auxiliary  equations  for  the  ether  components: 

(2.M  (f(z)     -     j(k^sin  a)M  [lp(z),       J(z)     -     j(-k^sin  a)/c4i^  V^(z) . 

The  quantity  L.(z)  may  be  thought  of  as  a  distributed  inductance  and  C.(z)  as 
a  distributed  capacitance,  whence  it  is  natural  to  define  a  characteristic 
or  wave  impedance. 


(2.?)         ^.(z)   - 


and  a  wave  velocity 


1/2 
L^(z)/C^(z)h    ,       i  -  P,n 


-1/2 
(2o6)         v^(z)  -  iL^(z)C^(z)^    ,       i  -  p,n. 

Corresponding  to  the  wave  velocity  of  (2.6)  we  have  a  wave  number 


-  u  - 


(2.7) 


k^(z)  =  oo/v^Cz), 


i  «  p,n. 


Using  these  definitions  we  have  in  place  of  (2.3) 


(2.6) 


dV^/dz  -  -j  Sk^(z)^^(z) 


r  I. 


di^/dz  -  -j  i 


k^(z)^^(z)f  V^, 


i  -  p,n, 


Since  the  quantities  V.  and  I.  are  solutions  of  a  pair  of  linear 
differential  equations,  they  may  be  regarded  as  the  components  of  a  two- 
dimensional  vector  (V, ,1. )o  In  order  to  describe  the  behavior  of  traveling 
waves  it  is  convenient  to  use  a  representation  of  the  vector  (V.,I. )  different 
from  that  implied  by  the  transmission  line  equations  (2.3).  The  new  vector 
representstion  will  be  connected  to  the  eld  one  by  a  linear  transformation: 


(2.9) 


V^(z)  -  U^(z)  +  D^(z) 


Ij^(2)  -  iUj,(z)  -  D^(z)j'  /  K^(z),      i  -  P,r, 


where  the  quantities  U.(z)  and  D. (z)  are  the  components  of  the  transformed 
vector  (U.,r. ).     The  transmission  line  equations   (2.8)  become: 


(2.1C)         dU^/dz     ♦  Sjk^  -  (d^^/dz)/2l^^^Uj  -  V(cS:^/dz)/2^^^Dj, 


dD^/dz     + 


-jk.    -  (d^./dz)/2>:^p^  -i-(d^^/dz)/2^.[u.,         i  -  p,n, 


If  we  assume  that  the  quantity  dJ^./dz  is  small  (in  some  sense  to  be 
specified  more  precisely  later  on),  we  should  expect  from  the  form  of  (2.30) 

that  U.  represents  a  wave  traveling  upward  only  while  D.  represents  a  wave 

traveling  downward  only.  It  is  clear  then  that  the  quantities  U.  and  D.  are  more 

appropriate  for  describing  traveling  waves  than  the  original  quantities  V. 

and  I . •  I 

1 
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3o   Approximate  solution 

In  view  of  the  assumption  that  di^./dz  is  to  be  small  it  is  natural  to 
attempt  to  solve  equations  (2.10)  by  means  of  a  perturbation  procedure.  If 
the  right  side  of  equations  (2.10)  were  to  vanish,  a  solution  for  z  <  A, 
except  for  two  arbitrary  constant  factors,  would  be; 

/   A 
(3.1)        U^(z)  -  ^p(z)expJjj  k^(s)ds 


D^iz)     '    ^p(z)  expJ-J 


z 
A 

k  (s)ds  V  ,         i  =  p,n. 


This  fact  suggests  a  change  of  variables  from  the  vector  (U.,D. )  to  the 
perturbation  vector  (u, ,d. ): 

(3.2)  U^(z)     -     u^(z)j:J/^(z)  exp^j   J   k^(s)ds 

Ti^(z)     -     d^{z)KY^{z)  exp  J-j       k^(s)ds    [>  ,  i  -  p,n. 

Equations  (2.10)  then  become: 

(3.3)  u^(z)  -  ^<^(z)/aCi(z)[^exp^-2j  I  k^(s)ds[.  d^(z) 

.A 


d^(z) 


-jr..'(z)/2J:.(z)l  exp  J  2J   k.(s)dsl  u.(2),    i  =  p,n. 


where  primes  denote  differentiation  as  usual,  Equations  (3.3)  have  solutions 
which  are  constant  if  the  medium  is  homogeneous. 

In  order  to  approximate  the  solution  to  (3.3)  hy  aeans  of  a  perturbation 
procedure,  it  is  helpful  to  change  (3.3)  to  an  integral  equation  by  integrating 
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both  sides.  We  have: 


(3.I4)    u^(z)  -  u^(A)  +  j  ^K[{y)/^^i7)\  exp  J-2J  J  k^(s)dsl  d^(y)dy 

d^(z)  -  <xl{7)/2K^{7)>  exp  ^2J   k^(s)dsl  u^(y)dy,        i  -  p,n, 


where  we  have  used  the  fact  that  only  the  upgoinp  wave  emerges  from  the  top 
of  the  inhomogeneous  layer,  and  therefore  d.  (/\)   «  0, 

We  may  obtain  an  approximate  solution  to  the  integral  equations  (3.U)  by 
a  standard  iteration  procedure.  A  single  iteration  can  be  obtained  by  sub- 
stituting the  inhomogeneous  term  u.  (/^)  and  the  term  d.  ( A)  which  is  zero  in 
(3.U)  into  the  integrals  in  place  of  the  solution  functions  u. (y)  and  d. (y). 
We  obtain,  after  replacing  u.(y)  by  u.(A)  and  d.  (y)  by  zero: 

(3.5)    u.(z)  -  u^(A) 

.A 


-3^(2)  ^   u^(A)   ^^i(y)/2^i(y)>exp  J2J  j  k^(s)d8ldy, 


1  -  p,n» 


h»       Iteration  of  a  general  system  of  integral  equations 

Consider  a  general  system  of  integral  equations  of  which  (3.I4)  is  a 
special  case; 

(U.l)    w,(z)  -  w,  (z)  +    7~  A.  .(y)w,(y)dy,  i-l,...,D. 


*  Hi 


1        i       ^T  iJ    J 


Let  w. (z)  be  the  m-th  approximation  of  an  iteration  procedure  for  solving 
equations  (Uol).  The  iteration  procedure  is  described  by  the  recursion  , 
relation: 
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(U.2) 


in 

W,(2) 


o 


f^n 


m-l 
2I  A. .(y)     w  (y)dy, 
j-1  ^ 


•j' 


i  *  1> •••)!)• 


The  solution  of  (I4.I)  is  given  by  the  limit 


w.(z)     -     lira      w  (z), 
itw>oo 


as  one  can  show,  using  a  standard  method.  The  error  coimnitted  by  stopping 
at  a  finite  value  of  m  can  be  estimated  at  the  same  time,  as  the  following 
demonstration  shows. 

From  the  recursion  formula  (h,2)  we  obtain  a  recursion  formula  for  the 
difference  between  two  successive  terms  of  the  sequence: 

(U.3)    w^(z)  -"wJCz)  -  j  f:   A.  j(y)  J  '"w^(y)  -"'w^(y)Uy,    i  -  i„ 


•  •» 


n. 


In  order  to  estimate  the  error  we  shall  need  a  bound  on  the  effect  of  the 
functions  A.  .(y)  in  (U.3).  We  therefore  define: 


n 


(U.M       ^(y)    -    ^X: 


III    ^j 


2  1/2 


We  must  obtain  a  single  estimate  for  the  largest  error  among  the  w^. 
From  ()4.3)  we  have 

1?       In 


<-'  &l"^'^'-"^^"lT-&)t'-<4" 


m-l           i"-2 
Wj(y)-    w^(y) 

dy 

xl 

2     7 

1- 

m 


-1  m-2 


hM  "j^y^"   ^J^^^ 


dy 


A 


n       n 


j»l  1^ 


)| 


w.(y)-  w  (y)  dy 

J  J  i 


^, 

'      n 

< 

iz 

/ 

V 

E  lAi^(y) 

i-l     ^-^ 


2 1  1/2  U     Im-l,    ,  m-2^    ^ 
5 J   "d^y^-  "j^y^ 


1/2 


dy 


/A 

.    .      iL      Im-l  m-2 


1/2 


dy. 
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With  inequality  (U.5)  we  can  obtain  an  estimate  of  the  difference  between 
two  successive  terms  of  the  iteration  sequence.  We  have  by  (U.^)  and  by 
induction: 


(U.6) 


& 


w^(z)-  w^(z)|  >   <  W 


A 


(' 


^'^^J^m       ^^(ym-l>^m-l-*-   ^^Vi^^Vv 


m 


^2 


where 


I    o 
W  -  max      -<7"  w.  (y) 

z  <  y  <  A  i^l ' 


1/2 


.  Through  successive  integrations  by 


parts  the  right  side  of  (h.6)  can  be  put  into  a  simpler  form.  Setting 

7 


S(y) 


M.(s)ds, 


we  have 


/A 


^A 


(U.7)    W  I  ^(y^)dy^  I  ^(y„.i)dy^.i  —  J  ^(yi)dyi 

.A 


^j  ^^^ym-i^^^Vi^^ym-i  j  ^'^\-2^^^^-2  -  j  ^^(yl)^l 

(W/2)J  Jd[s2(y^.^)]/dy„.J  dy^.^  j  ^(y,.^)6v,,2"j   ^^^^l^^^l 


A 

-  (W/3!)Ud 


^m-l 

A 


/      ym-2  ^2 


(W/iBl)  S"(A)  -  (W/ml)^   ^(s)dsV  . 


We  can  now  demonstrate  that  the  sequence  v.(z;   actually  converges. 


We 


have 
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N  0 


N-1 


m+1.    .        ni   /    \ 


Thus, 
(h.8) 


w^(z)   -  w^(z) 


lim 

N->oo 


N   (    ^ 
w^Cz) 


w^(z) 


m+1.    . 
w^(z) 


\i-)\. 


providing  the  limit  exists,  of  course.  We  have 


v^(z)  -  w^(z) 


w^(z)  -  w^(z) 


1/2 


Thus,  by  (U.6)  and  (U.7),  since 

(U.9)     Y.     (W/ni)  ^  i  t^(s)ds[>  -W^exp 
m»o 


A 
M.(s)ds 


1-  z 


n 


exists,  the  limit  in  (U.8)  exists,  and  the  sequence  w,(y)  actually  converges 
uniformly  in  the  interval  z  <  y  <  A,*  From  this  fact  and  (U.2)  it  follows 

that 

m        00,  .       .   . 
lim   w^(z)  =•  w^(z)  «  w^(z). 


m->oo 


where  w. (z)  is  a  solution  of  (U.l),  There  is  no  point  in  considerinp  the 
question  of  uniqueness  here. 

We  are  interested  most  in  the  error  committed  by  stopping  tho  sequence 
at  a  finite  value  of  m.  It  is  obvious  from  the  preceding  paragraph  that 
this  error  is  less  than  the  error  which  would  occur  in  (U.9)  if  the  infinite 
series  on  the  left  were  terminated  at  the  given  finite  value  of  m.  We  have 


for  the  error  e; 


(b.lO)    e,  -  w,(z)  -  w,(z) 


00 


m+j+1     m+j 


,  -,.,    .,..,. ^j    V-.      . 


z)-  w. (z)>,     i-l,..,,n. 


By  (U.IO)  and  a  repeated  application  of  the  triangle  inequality  we  have 
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(U«ll)       6     «     max  c. 

""       KKn    '     ^ 


<  ST^    w^Cz)-  w^(z) 


i=o     i=l  ^  1 


1/2 


We  combine  (U.6),   (U.7),   (U.9)   and  (ii.ll)   with  the  result  that 


(U.12)       e     <W^exp 


n(s)ds 

-2 


m 


|j.(s)ds 


<-  2 


7vi 


which  is  our  estimate  of  the  error  committed  when  the  sequence  is  stopped 
at  the  m-th  term. 


^,   Application  to  the  error  estimate  for  the  electromagnetic  problem 

We  are  interested  in  a  special  case  of  the  general  system  treated  in 
Section  h,   namely  the  integral  equations  (3.U)  associated  with  the  electro- 
magnetic problem  for  plane  stratified  media  discussed  in  Section  2,  We  can 
thus  identify  the  important  quantities  of  Sections  2  and  3  with  the  various 
quantities  of  Section  U  in  this  case: 

w^(z)  »  u^(2),      WgCz)  «  d^(z),    w^(z)  -  u^(A),  v^iz), 


A 
k^(s)ds)>  , 


A^^(z)  -  A22(z)  -  0,  A^2  '  <K^(^)/^K^Ml  ^^  i'^j 


Aj^^Cz)  -  ^^^(z)/25^^(z)j>  exp  J   2j   kj^(s)ds 


Thus,  by  (U.h)  and  (U,6) 

(5.1)         n(z)  -  (1//?)  k^(z)/^Jz)l, 
W  -  |u.(A)|   and   in  -  1. 

Thus,  by  (U,12),  when  (3.?)  is  our  approximate  solution  the  error  is  less  than 
e,  where 
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(?.2)    e  < 


u^(A) 


3Xp 


(1//?)   k^(y)/^.(y)|dy 


-1-(1//?)   k^(y)/^.(y)|dy 


The  integral  in  (^,2)   can  be  evaluated  more  explicitly.  We  have 


q  -  Iql  exp(j  |Z!j_), 


whera  0.  (y)  is  the  phase  of  S^^(y) .  Thus, 


and  therefore 


l</^ili  ikilV  Iqll  *  l<l- 


We  therefore  have 

/A 

(5.3) 

ki(y)Ai(y)ldy  <  1    llqlV  Iq  l|dy  *  J    I0i|dy. 


Now  let  the  numbers  z     be  all  the  points  at  which    |^. |     changes  sign  in  the 
range  of  integration,  except  for  two  values  of  v  at  which  2     is  equal  to  z 
and  A,*  and  let   f   be  the  larger  of  the   two  ratios 

We  define  a  quantity  5^.   which  might  be  called  the  total  variation  of  >^^,   to 
be  the  product  of  all  the    V    : 


5^. 


TTL. 


Also,  let  the  numbers  z,  be  all  the  points  at  which  <j/.  changes  sign  in  the 
range  of  integration  except  for  two  at  which  z,  equals  z  and  z2l»  Then,  by 
(5.3), 
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(^<.U)       I  k{(y)/^i(y)|ay  <  log  (eq)  *  r  |9f(z^_^^)  -  0(z^)  |. 


The  expression  (5.2)  can  now  be  written 

,(1//?) 


(5.5)     e<  |u^(A)K(6^.)^-^/  ''"^  exp 


(1/y?)  y  |9r(z^^^)-0.(2^)| 


i+(i//?)iog(6^)+(i/v^)  r;  Mi(2^^i)-  0i(z^) 


,1 

The  case  of  zero  conductivity  is  the  case  in  which  0.    is  zero.  Then  we 
have  simply 

(5.6)        c  <  |u^(A)|  ^(5^^)      -l-(l/-/?)log(6^)y  . 


Ihe  inequality  (5.6)  can  be  replaced  by  a  cruder  but  simpler  estimate 
of  the  error.  We  have  for  x  >  0 

1/  (x+1)  >  l-xj 

after  integrating  both  sides  we  have 

(5.7)        log  (x+1)  >  x-x^/2. 
From  ('^.6)  and  (5.7)  we  obtain 

(5.6) 


<  \n^(A)\}itK)^^^^^'lV/2, 


The  exponent  1/  ■/?  in  (5.6)  can  be  left  out  in  a  slightly  cruder  estimate. 

6.   Frror  estimate  for  the  reflection  coefficient 

We  are  more  interested  in  computing  the  reflection  coefficient,  which 
is  the  ratio  of  D.  to  U.,  than  we  are  in  the  actual  field  quantities  as 
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determined  by  an  initial  value  problem.  This  means  that  we  should  also  in- 
vestipate  the  error  in  the  reflection  coefficient  due  to  the  existence  of 
an  error  in  the  field  quantities.  The  reflection  coefficient  R.  is  given  by 

(6.1)  R^  -  D^/a^. 

Let  D.  A^.   be  the  approximate  vailues  of  D.  and  U,  based  on  the  approximations 
(3.^)  for  d.  and  u, o  Let  e.  be  the  difference  u,-u.  and  6-  the  difference 
d.-d,  o  Let  C.  /^,  be  our  approximate  reflection  coefficient.  We  have  for 
the  error  in  computing  R . ,  assuming  zero  conductivity, 

(6.2)  .    iVio-^i^iol/IViol   •    l(^io*^2^^o-(^io*^l)^iol/IViol 

-    1^0^2-^io^ll/IViol   -    1^2/^i-^io  V"il 

<    le^l/luJHRiollhl/l^l  ^   l^2l/|l^ioH«ll 

Since  every  term  in  (6.2)  contains  or  is  bounded  by  a  quantity  containing 
u.  (^)  as  a  factor  we  can  set  u.(^)  ■  1,  Then 


(6.3)  Cj.  <   Icgl    /  , 


l-|e^|U    iRiJIeJ   /|l-|ej|  . 


A  difficulty  seems  apparent  in  (6.3)  in  case  the  approximate  reflection 
coefficient  R.  has  a  magnitude  of  the  same  order  as  the  error  62 »  for  then 
the  error  term  may  be  large  enough  to  obscure  any  significant  approximation 
of  the  reflection  coefficient.  It  is  therefore  necessary  to  determine 
whether  | e- |  is  sufficiently  small  to  justify  assuming  that  R.  is  zero  in 
this  case. 

We  observe  that  in  the  electromagnetic  problem  the  matrix  A.  .(z)  of  our 


system  of  integral  equations  has  only  zero  elements  in  its  main  diagonal. 
The  product  of  tvro  such  matrices  is  a  diagonal  matrix.  It  follows  that  if 
we  perform  another  iteration  of  (3.h)  to  obtain  a  correction  tc  (3.5),  the 
approximation  to  the  component  d. (z)  remains  unchanged  while  u. (z)  gets  a 
correction  term  added  to  it.  Our  conclusion  is  that  the  estimate  of  the 
error  in  our  approximation  to  d. (z)  is  of  hipher  order  than  we  have  actually 
assumed,  corresponding  in  fact  to  one  more  iteration  of  the  integral  equation, 
although  our  error  estimate  for  the  component  u.(z)  is  more  or  less  correct. 
Thus,  jcj,!  is  bounded  by  a  quantity  of  higher  order  than  the  bound  on  |6-|o 
Thus,  according  to  (6.3),  the  error  e  will  either  be  of  higher  order  than 
R.  or  else  we  may  set  R.  equal  to  zero  with  an  error  of  higher  order  than 
|e-|.  Since  |e^  |  is  of  hiphpr  order  than  the  aniplitude  of  the  incident  wave, 
we  conclude  that  our  approximation  to  R.  is  significant,  and  a  significant 
bound  on  the  error  committed  in  this  approximation  is  given  by 

(6.U)        e^  <  e^/2(l-e)  +  |R^^|e/(l-e) 

where  e  is  given  by  (5. "7)  or  (5.6),  Implicit  in  this  result  is  the  assumption 

that  c  is  less  than  one  as,  in  fact,  it  must  be  in  order  that  the  approximation 

have  any  real  meaning. 

When  the  conductivity  is  different  from  zero  the  error  estimate  (6.U) 

should  be  multiplied  by  the  exponential  factor 

A 


(6.5) 


exp  <-j   k.(s)ds 


1 
z 


which  is  obviously  one  when  k.(s)  is  real  in  the  interval  from  z  to  ^, 

7,   Summary  and  conclusions 

We  have  treated  approximately  the  reflection  problem  arising  from  a  plane 
electromagnetic  wave  incident-  on  a  plane  stratified  inhomogeneous  aiedium,  and 
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we  have  obtained  a  bcnjnd  on  the  error  committed  in  the  approximation.  The 

perturbation  method  employed  in  this  treatment  can  be  characterized  roughly 

by  saying  that  it  holds  for  a  medium  of  small  variation,  A  precise  meaning 

is  given  to  the  term  "  small  variation"  if  we  define  it  to  be  equivalent 

to  the  statement  that  the  quantity  5^.  (defined  in  Section  5  to  be  the 

total  variation  of  Y.  )   i5  small. 

1 

Two  rather  unusual  features  of  the  approximation  treated  here  should 
be  noted.  First,  it  is  valid  for  all  frequencies  co  to  within  the  error 
e  (defined  in  Section  6)  which  is  independent  of  frequency.  Second,  it  is 
valid  even  when  derivatives  of  the  electrical  parameters  e,  \i.   and  cr  of  the 
medium  have  jump  discontinuities. 

For  easy  reference  a  summary  of  our  treatment,  including  the  relevant 

formulas,  follows: 

/I 

( 1 )  We  assume  a  cartesian  coordinate  system  whose  z-axis  is 

normal  to  the  planes  on  which  the  electrical  parameters  are  constant. 
The  x,y^lane  boxinds  the  inhomogeneous  medium  below,  while  a  plane 
z  ■  ^  bounds  it  above.  The  y,z-plane  is  the  plane  of  incidence,  and 
the  angle  of  incidence  is  designated  by  g  (see  Figure  1). 

(2)  Letting  the  subscript  p  stand  for  polarization  parallel  to 
the  plane  of  incidence  and  n  for  polarization  normal  to  the  plane  of 
incidence,  we  write  for  the  field  vectors  E  and  H 

E-E+E,  H-H+H, 

p    n*  p    n' 

where  the  vectors  have  the  form 

E   -  (0,V  ,  £)  exp(-jk^y  sin  a), 
H   -  (I  ,0,0)  exp(-jk^y  sin  a), 
^n  ■  ^\'°»°^  exp(-Jk^y  sin  a). 


"n  "  ^°»^n>^^  exp(-jk^  sin  a). 
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with  k     -  CO  yTir  ,     Inductances  L  ,  L^  and  capacitances  C^,  C^  are 


defined  by 


L  (z)  -  (u  e^sin  a)/  t{z)-\i{z), 
p  o  o 

L  (z)  »  |J.(z), 
n 

Cp(z)  -  --feCz), 

C   (z)  -  e(z)   -  (u  c^sin^a)/^(z). 


n 


o  o 


(3)     Fouations   for  the  field  vector  components  are 
dV^/dz     -     -jaiL^(z)l^, 
dl^dz     -     -jc^^(z)V^,  i  -  P,n, 

along  with 

6(z)     -     nV^  a)/co^  ^    yz). 


^(z)     -    1(-Vin  a)/co|^  r  V^(z). 

(U)     The   characteristic  irepedances  K^i^)  are  defined  by 

1/2 


q(z)     -    '^L.(z)/C.(z) 


P>n, 


and  the  wave  numbers  k^(z)  by 


k^(z)     -     CO  JL^(z)C^(z)p/^  i  -  p,n. 


(5)  The  upRoing  components  Uj^(z)  and  the  downgoijig  components  D^(z) 
of  the  traveling  wave  representation  are  determined  by  the  relations 
V^(z)  -  U^(z)  +  D^(z), 
I^(z)  .  ju^Cz)  -  Di(z)j.  /q(z),         i  -  P,n. 
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(6)   The  perturbation  factors  u.(z)  and  d.(z)  for  these  quantities 
are  given  by  the  relations 

U^(z)  -  u^(z)  ^'•^^(z)  expj  j  j  k^(s)d! 

^   z 

D^(z)  -  d^(z)  ^^/^(z)  expj  -j  J  k^(s)ds|  ,       i  -  p,n. 


(7)  The  approximate  solution  is  given  by 
u^(z)  ^  u^(A) 

d^(z)  ^    u^(A)  jx:^i(y)/2^i(y)exppj  j  k^(p)dslUy,   i  -   p,n. 


(8 )  Let  the  numbers  z  be  the  points  at  which  the  derivative 
|s^.  (z)l'  changes  sign,  except  for  two  values  of  z  which  are  z  and 
ZA  and  let  P  be  the  greater  of  the  two  ratios 

Mvi^i  /  l^i^V  ^^    iq(s)i  /  tq(vi)i- 

Then  the  total  variation  6^.  of  ^.  is  defined  to  be 

^i   •   TT  Tv . 


(9)  Let  0.  (z)  be  the  phase  of  5^.  (z)  and  the  numbers  z,  the  points 
where  the  derivative  0. (z)  changes  sign,  except  for  two  values  of  z, 
which  are  z  and  /\,    A  bound  on  the  error  e,  in  our  approximation  to 
the  vector  (u. ,d . )  is  given  by 


/ 


< 


1^(1/  y?)(log  6^^  +  rii^i^'^x*!^'  "^i^^X^I^I*      i  -  P,n. 
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(10)  If  we  assume  that  the  conductivity  is  zero  then 

<  (5^^)^^/  /^^  -  Jl  +  (1/  /?)log(6q)i  ,       i  -  p,n. 


"i 


a  cruder  estimate  is 


MUi)<'/'^'4%:k4'/^. 


i  -  p,n. 


(11)     A  bound  on  the  error  estimate   e     for  the  reflection  coeffi- 

r — 

cient  is 

Bj.^  <^e^/?(l-e^)+|R^^|e^/(l-6^)Vexp  J-j    1  k^(s)ds>,       i  »  F,n, 

where  R,     is  our  approximate  value  for  the  reflection  coefficient 


determined  by  the  ratio 


So  '  ^io/^io»  i  "  P»" 


of  oxir  approximations  D.  and  U.  to  D.  and  U . , 
'^^ 10  —  10  —  1  1 
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